Numerical Analysis of Boosting Scheme for Scalable NMR Quantum Computation 
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Among initialization schemes for ensemble quantum computation beginning at thermal equilib- 
rium, the scheme proposed by Schulman and Vazirani [L. J. Schulman and U. V. Vazirani, in 
Proceedings of the 31st ACM Symposium on Theory of Computing (STOC'99) (ACM Press, New 
York, 1999), pp. 322-329] is known for the simple quantum circuit to redistribute the biases (po- 
larizations) of qubits and small time complexity. However, our numerical simulation shows that the 
number of qubits initialized by the scheme is rather smaller than expected from the von Neumann 
entropy because of an increase in the sum of the binary entropies of individual qubits, which indi- 
cates a growth in the total classical correlation. This result — namely, that there is such a significant 
growth in the total binary entropy — disagrees with that of their analysis. 

PACS numbers: 03.67.Lx, 05.30.Ch, 05.20.Gg 



I. INTRODUCTION 

Initialization is an indispensable stage in ensemble 
quantum computation beginning at thermal equilibrium, 
especially in NMR computations [E 0, 0, S H • For a 
model of the NMR system, it is customary to consider 
a mixture of molecules with n spins- 1/2 under a strong 
static magnetic field as shown in Fig.^on the assumption 
of the independence among initial spins; local couplings 
among spins in a molecule are negligible compared with 
their Zeeman energies Ea- The original bias (polariza- 
tion) e of a spin at temperature T is given by 
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where ks is Boltzmann's constant. This spin is in ||) 
with probability (1 + e)/2 and in | J.) with probability 
(1 — e)/2. We regard the up spin as the bit of and the 
down spin as the bit of 1. Quantum operations are per- 
formed in this n-qubit 0, Q ensemble system. However, 
it contains 2" component states |0 • • ■ 0), . . . , |1 ■ • ■ 1), un- 
less T = 0; we have to extract the final state evolved from 
a particular component state to obtain the outcome of a 
quantum computation for any standard quantum algo- 
rithm. Suppose that the correct input is only the state 
|0}® n . While its original population is largest among all 
the component states, the product of probabilities of each 
bit being is small for large n. The more bits we need, 
the larger biases are required to extract the final state. 
As it is impossible to achieve £ ~ 1 by cooling sample 
materials, algorithmic initializations 0, 0> 0, S Ell EI 
are utilized to prepare an initialized state: namely, a 
particular component state with an enhanced signal. 

Standard quantum algorithms are designed to work 
with a pure-state input. In contrast, there are a few 
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FIG. 1: A model of the NMR system for spin-1/2 nuclei. 



quantum algorithms for the model of deterministic quan- 
tum computation with a single pure qubit (DQC1) [T^ 
and similar models, including one for prime factorization 
[l3| (an extension of Shor's algorithm 0]); those are de- 
signed to work with a highly mixed state, even with an 
input string that comprises a single pure qubit (namely, 
an initialized qubit) and the rest qubits in a maximally 
mixed state. As the family of DQC1 models does not 
provide all the replacements of standard quantum algo- 
rithms, initialization of many qubits for ensemble quan- 
tum computers is still of importance. 

There are two types of initialization methods: one is 
averaging (e.g., effective pure states @); another is data 
compression. The latter has an advantage in respect of 
scalability because the former needs resources growing 
exponentially in n while the latter does not. The boost- 
ing scheme proposed by Schulman and Vazirani 0, ll(| is 
one of the latter type of methods. Actually, their scheme 
is a sort of polarization transfer which has been com- 
monly used and thoroughly studied in the realm of NMR 
spectroscopy as pointed out in Ref . ^jjj . There have been 
notable studies of the general upper limit on a polariza- 
tion enhancement in a polarization transfer (e.g., Refs. 
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[l6l IT?! Hsf ) . The principle of the boosting scheme is the 
bias-enhancing permutation as we will see in Sec. |H] In 
the scheme, a quantum circuit @ properly composed of 
basic bias-boosting circuits is used for redistributing the 
biases to generate the block of / cold qubits with biases 
greater than e co id such that 

Prob{0 1 2 • • • 0*- X ; } > (l±^±\ ' (2) 

is large enough to distinguish 10)®' among the signals of 
2 l component states. Their analyses 0, El resulted in 
that a compression down to almost the entropic limit is 
possible (in the limit of large n) in 0(n log n) time in the 
scheme to generate n(l — o(l)) — S initialized qubits each 
of which has bias 1 — exp(— n 6 ^ 1 )) (which converges to 
1 as n becomes large), where S is the von Neumann en- 
tropy 0] of the ensemble However, this is incorrect 
unless the original biases are in the neighborhood of or 
1 Irrespectively of n, the condition of original biases 
< e 1 leads to the circumstance that the probability 
of a spin being up has classical correlations (classical de- 
pendences) with those of other spins within several steps 
after the beginning of the scheme. Owing to the correla- 
tions, I can be much smaller than n — S. Especially for 
the practical setting of n < 10 3 , an accurate evaluation 
of the scheme is of importance since a large mean bias 
such as (e) — 0.4 0] is required to obtain a sufficient 
number of initialized qubits 20]. 

In this Paper, we report a numerical analysis of the 
boosting scheme based on a molecular simulation. Sev- 
eral features of the scheme relating to the total binary 
entropy have been investigated by using the simulation 
with the aid of matrix calculations. It has revealed the 
fact that the number of qubits initialized by the scheme 
is rather smaller than expected from the von Neumann 
entropy even for large n. This is not a surprise; it is quite 
possible according to the theory of macroscopic entropy 
described in Sec. II I II 

II. BOOSTING SCHEME 

Enhancing the bias £ 2 ; of the ith spin is an opera- 
tion equivalent to increasing the population of a set of 
molecules whose ith spin is up [si = £i(t) where t de- 
notes time step] . This is a general conception of initial- 
izations for mixed states. As we consider a large number 
of molecules under a strong static magnetic field such as 
those in a common liquid- or solid-state NMR system, an 
n-qubit thermal equilibrium state can be written as the 
density matrix 

Pe q =p7®p7®---®p e n U®p e n q , (3) 

where p^ q is the density matrix representing the original 
state of the ith spin with the bias £i(0): 

*■-[ I 4m)- w 



Let Cfc = Cfe (t) denote the population of component state 
\k) at time step t. Then it is also represented as 

2"-l 
fc=0 

= diag{c & (0),ci(0),... ) ca»_a(0),C3»-i(0)}, (5) 

where Cfe(0) are the original populations. The bias of 
the ith spin is enhanced by any unitary operation that 
permutes a component state with a large population to 
some component state whose ith bit is 0, since such an 
operation increases the population of a set of component 
states whose ith bit is 0: 

Pi= Ck - ( 6 ) 

A: | ith bit is 

This population is identical to the probability that the 
ith bit is 0. Because a component state which consists of 
a larger number of 0's and a smaller number of l's has a 
larger population among the 2™ component states in most 
cases, any unitary transformation that maps a sufficient 
number of component states with more than n/2 0's to 
those with 0's in specific I' bits enhances the biases of the 
I' spins. The initialization of I > I' qubits can be realized 
by a combination of such transformations. When the 
transformation U is only a product of permutations of 
component states, as it is nothing but the exchanges of 
values of Co, ... , opt-x, the evolved density matrix p — 
UpeqU^ is still represented as 

2 n -l 

p= £>|fc)<fc|, (7) 

fc=0 

although it is often impossible to write p in the form of 
Eq. ©. 

In Schulman and Vazirani's boosting scheme (or sim- 
ply called the boosting scheme), the total operation for 
an initialization can be composed of the 3-qubit basic cir- 
cuits shown in Fig. ^3] or the 4-qubit circuit shown in 
Fig. which perform permutations of component states 
as mentioned above. It should be noted that these cir- 
cuits were not distinctly written in the original descrip- 
tion of the boosting scheme. Either the 3-qubit circuits 
or the 4-qubit circuit can be a basic operation of the 
scheme. If one interprets the original description as it is, 
one may have the 4-qubit circuit. We, however, use the 
3-qubit circuits rather than the 4-qubit circuit because 
the 3-qubit circuits are better with respect to the rate of 
initialization. We do not write details about it here, but 
we will describe a reason in Sec. IIVI From now on, we 
will regard only the 3-qubit circuits as the basic circuits 
of the scheme. 

The operation of the basic circuits is illustrated in the 
truth table of Table[Q By the circuit of Fig.^a), a com- 
ponent state with a large population — i.e., a component 
state in which two of the three qubits are — is mapped 
to a component state in which the first of the three qubits 
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FIG. 2: Basic quantum circuits for constructing an initializa- 
tion circuit in the boosting scheme. The CNOT gate performs 
a bit flip (namely, a NOT operation) on the qubit c when 
b — 1. The Fredkin gate swaps the qubit a with the qubit b 
when c — 1. 
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FIG. 3: Another boosting circuit with a 4-qubit input. The 
CC-SWAP gate swaps a with c when bd — 11. 



is 0. The probability that the first qubit a is after this 
basic operation is 

Prob{a = 0} out = Prob{000} m + Prob{001} m 

+Prob{010} in + Prob{100}„, (8) 

where, for given qubits a, 6, and c (from upper to lower 
in the circuit) with their binary values x a , Xb, and x c , 



Prob{x a a; b a; c } = c k . 

k\abc—x a x^x c 



(9) 



Suppose that no correlation exists among P a ,Pb, and P c 
and that the biases e a , £&, and e c are positive, in the 
input. Then the bias of the qubit a in the output of the 
circuit is 



E 'a = 2 + £ fc + £ c — £ a £fc£c) ; 

and e' a becomes larger than e a under the condition: 

% + £ c 



1 + SbSc ' 



(10) 



(11) 



Similar calculations give the biases for b and c in the same 
output: 



Sb — e c + e a Sb£ c ) ■ 

- E-b^-c- 



(12) 
(13) 



After the circuit of Fig.^a), the bias of the qubit b may 
be smaller than 0; the circuit of Fig. ^b) is added to 
invert the bias in this case. In the case where three spins 



TABLE I: Truth table for the basic circuits in Fig. |5] for the 
boosting scheme. Here, the input of (b) is the output of (a). 

3-qubit boosting 



Input 


Output of (a) 


Output of (b) 


000 


000 


010 


001 


001 


011 


010 


011 


001 


011 


100 


110 


100 


010 


000 


101 


101 


111 


110 


111 


101 


111 


110 


100 



are polarized uniformly with the bias e un i in the input, 
Eq. JUII) leads to 

e a = 2 £uni ~ T^-uni- (14) 

Thus, the basic circuits can be used for boosting the 
biases of target qubits. It is expected that one can 
achieve the component state |0)®' with a very large pop- 
ulation simply by collecting initialized qubits after suffi- 
cient steps of application of these circuits to every sub- 
group which consists of three qubits that satisfy inequal- 
ity Ijlljl (although the initialization rate is another sub- 
ject). Equations (fT77f> . fl2} , and (fT^f may be used for 
a rough estimation of the behavior of a quantum circuit 
composed of the basic circuits. 

However, the designing of a whole circuit in the boost- 
ing scheme is not so simple. At every step in the scheme, 
one recreates subgroups by picking up every three qubits 
from the n qubits in the order of their biases, from larger 
to smaller, so that the condition of inequality pi|l is sat- 
isfied in almost every subgroup. It is better to avoid 
applying the basic circuits to those in which the condi- 
tion is not satisfied or undo it afterward. In order to 
recreate the subgroups, one has to forecast the biases of 
qubits correctly, which requires a more precise estimation 
owing to the following reason: After several steps in the 
scheme, the bias of a qubit is not independent of those 
of others, i.e.. 



Prob{ij = 00} ^ PiPj 



(15) 



for a pair of qubits in the n-qubit string. This is 

because the ensemble after a CNOT operation on 
is not determined by the biases (si,Ej) alone but depen- 
dent on the populations {cfc}. The exception is only the 
case where the original biases are very close to or 1. 
Therefore, the use of independent probabilities for repre- 
senting spins in constructing the circuit induces inaccu- 
racy for the large size and/or large depth of the circuit. 
Each molecule has to be dealt with for correct forecasts 
of biases. 

Moreover, even with the forecasts of biases, it is quite 
possible that the bias of the qubit a is not boosted in 
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FIG. 4: A unitary compression operation U for redistributing 
the biases to create the block of qubits with small binary 
entropies and that with large binary entropies. 



some subgroups at almost every step except for several 
beginning steps in the scheme. This is due to classical 
correlations among qubits. The condition of inequality 
is not an accurate condition for the success of bias 
boosting in a subgroup of correlated qubits. It is nothing 
but a rough estimation of a correct condition. Thus, 
one ought to undo the operations that fail in enhancing 
the biases of target qubits so as to achieve better bias 
enhancements. The details of the circuit designing based 
on the forecasts and undo operations are described in 
Sec. [V] and Appendix IU1 

Note that even the circuits that succeed in enhancing 
the biases alter the amount of correlation. Although the 
three qubits in each subgroup are nearly equally biased, 
this does not mean the binomial distribution, i.e., it is 
not relevant to the independence of the qubits; their bi- 
ases can be correlated [their individual probabilities given 
by Eq. ((SJ) can be dependent] to each other. The classi- 
cal correlations affect the efficiency of initialization — this 
"efficiency" is not the efficiency in running time nor that 
in working space — as described in the next section. 



e l = -1 + 2 



E 



Cfe. 



(17) 



felith bit is 



On the other hand, when there can be nonzero off- 
diagonal elements of the density matrix for the system, 
we have to consider the intrinsic bias defined below rather 
than the superficial bias that is calculated only from pop- 
ulations. The intrinsic bias of a qubit is a quantity ob- 
tained from a unitary diagonalization of the reduced den- 
sity operator of the qubit. The reduced density operator 
of the ith qubit can be given by the partial trace 



E 



(18) 



where p — p(t) is the density matrix for the whole en- 
semble; {IV'fc)} = {IV'o); • • • > IV^-i-i)} 1S a complete or- 
thonormal system of the state space of the rest qubits, 
such as 



{|0i0 a ■ • • 4 _iO<+i • • • w _i0„>, 
|0i0 2 ---0 i _i0 <+ i---0 n _il n ), 

3 

1 1 1 1 2 • • • li-llj+1 • • • lri-lO n ), 
1 1 1 1 2 ' ' ' li-llj+1 ' ' ' Iri-lln)}- 



(19) 



Let d p~i be the matrix with the diagonal elements Ci$ 
and Cf i obtained from a unitary diagonalization of pi 
such that Ci t o > c%,i > 0. Then the intrinsic bias of the 
ith qubit is defined as 



III. CLASSICAL CORRELATIONS HIDDEN IN 
MACROSCOPIC BINARY ENTROPIES 

The initialization of an NMR quantum computer with 
a method that uses a data compression is equivalent to re- 
distributing the binary entropies of qubits (the entropies 
of macroscopic bits), {Hi = ^4p- log 2 j^— + ^—^ L log 2 

t— — }, to create two different blocks of qubits: one has 
a very small entropy close to and the other has a very 
large entropy as illustrated in Fig. The block of I 
qubits with very small binary entropies is used as |0)® z . 
In this respect, an important entropy measure is the sum 
of the binary entropies of individual qubits: let the effec- 
tive entropy 

1=1 v 

(16) 

where £j = £i(t) is the bias of the ith qubit at time 
step t. When the state of the whole n-qubit ensemble 
is represented as a diagonal density matrix, it is given by 



Cifl — c%,\ 
d Ci,0 + d Ci,i ' 



(20) 



This definition of intrinsic bias is valid in the sense that 
its value is independent of a macroscopic spin direction. 
The direction of the ith macroscopic bit pi can be mod- 
ified even with only single-qubit rotation gates. If we 
restrict the subsequent quantum gates to single-qubit ro- 
tations at a particular time step, 1 ~^ Ei is the maximum 
achievable probability of the ith qubit being after that 
step. Thus the intrinsic bias is a plain measure of pu- 
rity of a qubit. Since an initialization is to produce pure 
qubits, we regard £$ as the actual bias of the ith qubit 
denoted by e% in the definition of S e : namely, in Eq. i|lbfl . 
When completely initialized qubits are required, S e is a 
lower bound for the length of the compressed string at a 
particular time step in an initialization process; n — S e 
is an upper bound for the number of initialized qubits at 
that step |45j| . S e varies step by step during the initial- 
ization. 

Now we will see an upper bound for the number of ini- 
tialized qubits under a realistic condition. Suppose that 
we do not need completely initialized qubits. Then ini- 
tialized qubits are allowed to have small binary entropies 
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(Hj) whose average value (y J2j=i Hj) ^ s equal to or less 
than some real number a (0 < a « 1). The total binary 
entropy of / initialized qubits is equal to or less than al. 
The block of rest n — I qubits (namely, the compressed 
string) therefore has the total binary entropy E > S e —al. 
Hence the following inequality holds: 

n-l>E>S e -al. (21) 

This leads to 

l<-^—(n-S e ). (22) 
1 — a 

Setting (3 — , we have 

l<(l+P)(n-S e )<n(l-S e /n + 0). (23) 

In an initialization, even if there is no need of com- 
pletely initialized qubits, one picks up I cold qubits at the 
terminal step of it in order to achieve Prob{0i • • • 0;} > c 
(0 <C c < 1 is a real number). The value of c is chosen 
by one's need. Accordingly, a limitation of the value of 
a can be estimated as follows: Suppose that the density 
matrix for a state after initialization has no nonzero 
off-diagonal element. Let us assume for a moment that 
one obtains I initialized qubits with independent biases 
and the rest qubits with biases that may be correlated. 
If the joint probability of the initialized qubits being 
zeros is Prob{0i • • • 0;} > b (0 < b < 1), then the 
binary entropies Hj for the qubits satisfy the inequality 

E l i=iHj < -^ 1/i log 2 & 1 /' + (i-&i/')log 2 (l -b 1 ' 1 )], 
which is easy to prove by mathematical induc- 
tion. From the contraposition, if Ylj=i Hj > 
-Z[6 1 /'log 2 6 1 / i + (l-6 1 /')l g 2 (l-6 1 /')] 5 then 
Prob{0i ■ ■ - 0;} < b. Therefore, on the present assump- 
tion that the biases of initialized qubits are independent, 
a necessary condition for achieving Prob{0i • • • 0/} > c is 
that 

a^-c^log^-Cl-c^log^l-c 1 /')- (24) 

Hence on the present assumption, a tends to as I in- 
creases, and so does /3. Possible values of a and (3 are 
estimated with inequality l|24l) by specifying a value of c 
and assuming that I is in a realistic range. For example, 
when one needs c = 0.99 and assumes I > 16, the values 
of a and are estimated to be at most 7.59 x 1CP 3 and 
7.64 x 10 -3 respectively. Indeed, it is, however, possible 
that a does not tend to with increasing I to achieve 
Prob{0i ■ • - 0/} > c when there are correlations among 
initialized qubits, but even for such qubits, of course the 
condition a < — clog 2 c — (1 — c)log 2 (l — c) is neces- 
sary. Moreover, in the case where one needs the bias 
of an initialized qubit to converge to 1 as n grows — e.g. 
£coid = 1 — exp(— n x ) (constant x > 0) — a and f3 tend 
to as n increases — i.e., (3 = o(l) in this case. As we 
have seen, inequality l|23|l gives an upper bound for the 
number of initialized qubits at a particular time step for 



a realistic initialization. In other words, S e — f3(n — S e ) is 
a lower bound for the length of the compressed string at 
that step. /3 is a small number, which is estimated to be 
less than 10 -2 in most of realistic initializations to make 
several bytes almost pure. 

Generally, it is believed to be possible to almost faith- 
fully compress the string of an n-qubit ensemble, repre- 
sented by density matrix a, down to the entropic limit 
with the rate S(a)/n for large n in all kinds of quantum 
computer including NMR computers (as long as n can be 
large) because of the analogy to the quantum noiseless 
coding 0, where 

S(a) = -Tr(alog 2( r) (25) 

is its von Neumann entropy. It is desirable, and be- 
lieved to be possible, to generate a quantum circuit to 
accomplish it within the depth of the circuit increasing 
polynomially in n. Indeed any unitary transformation 
of a density matrix conserves its von Neumann entropy, 
but it would be probably impossible for a compression 
scheme which generates a quantum circuit according only 
to macroscopic quantities of individual qubits, e.g., bi- 
ases, to do such an ideal compression for an arbitrary in- 
put state. For the initialization of NMR quantum com- 
puters, this is somewhat clear when we rewrite S e by 
using the reduced density operators {pi} of individual 
qubits. Regarding as the actual bias of the ith qubit, 
we have Hi = S(pi). Consequently, S e is rewritten as 

n 

S e = J2S(Pi). (26) 
i=i 

From this notation, it is clear that 

S e -S = S(p || pi® ■••<8>Pn), (27) 

where S = S(p), and S(p || pi ® • • • 8 p n ) is the Umegaki 
relative entropy [2l| of p with respect to pi ® ■ ■ • ® p n . 
In general, the Umegaki relative entropy between two 
density operators ci and cr 2 is defined as 

5( ( r 1 ||a 2 ) = Tr^ ( 7 1 log 2 ^^ . (28) 

It is a measure of the distance (or dissimilarity) between 
the two operators. On the right-hand side of Eq. 127fl . 
p represents the whole ensemble while p\ ® • • • (8 p n is a 
tensor product in which correlations are freed (i.e., cor- 
relations among pi , . . . , p n are obliterated in taking the 
tensor product). Hence S e — S is a measure of total cor- 
relation [2^ of p: namely, a measure of total correlation 
of each qubit to the other qubits. Obviously, can 
be much larger than S when there is a growth in the 
total correlation during the initialization. There is no 
known way to avoid the creation and growth of corre- 
lations among biases when a data-compression circuit is 
designed by using solely the values of biases at every step 
instead of those of commonly-used quantities such as the 
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(approximate) populations of possible bit sequences [23| . 
the (approximate) length and frequency of patterns in 
microscopic strings, etc. 

Let us consider the boosting scheme, in which every 
operation is a permutation of diagonal elements in a di- 
agonal density matrix which is originally p eq as we have 
seen in the previous section. Then, clearly, there is no 
entanglement during and after a bias boosting process; 
S e — S is a measure of total classical correlation hidden 
in the macroscopic string (the string which consists of 
n macroscopic bits p\,...,p n ). Thus, in the boosting 
scheme, an increase in S e measures a growth in the total 
classical correlation. 

In addition, the maximal efficiency of initialization of 
the macroscopic string is defined by 



where S^ nd is the value of S e after the application of an 
initialization scheme. This is the ratio of an upper bound 
for the number of completely initialized qubits that can 
be achieved by the scheme to the one that is expected 
from the von Neumann entropy. Similarly, the maximal 
efficiency of compression of the macroscopic string is de- 
fined by 

r c = S/St nd . (30) 

One can evaluate the performance of the scheme by cal- 
culating these efficiencies. 

IV. GROWTH IN THE AMOUNT OF 
CORRELATION FOR A BASIC BOOSTING STEP 

As we have seen in Sec. IIIII an increase in the amount 
of correlation must be suppressed to achieve a high ini- 
tialization rate. In this section, we therefore evaluate the 
amount of classical correlation generated by a single step 
of basic boosting circuits. 

Suppose that there are three qubits a, 6, and c with 
uniform bias e at thermal equilibrium. For this initial 
state, the effective entropy defined by Eq. (|16|l is equal 
to the von Neumann entropy S. Now we apply the 3- 
qubit circuit of Fig. EJa) to the qubits. The biases of 
the qubits become 36 ~ g , eH 2 £ , and e 2 , respectively. By 
this operation, classical correlations are induced, and the 
effective entropy of the output state, S° ut , is greater than 
S as shown in Fig. The discrepancy between S° ut 
and S is rather serious for the highly but not completely 
polarized qubits because S° ut — S is fairly large in despite 
of the small values of S. 

Similarly, we can evaluate the amount of classical cor- 
relation generated by the 4-qubit circuit shown in Fig. [3] 
When there are 4 input qubits a, b, c, and d with uni- 
form bias e at thermal equilibrium, the biases in the out- 
put of it are 3e ^ e , e 2 , and e 2 , respectively. Thus 
the effect of bias enhancement looks almost the same as 
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Original bias £ 

FIG. 5: Plots of ST*, S, and S° ut - S for 3 qubits as func- 
tions of e, the uniform original bias of the qubits. S° u is the 
effective entropy of the output state after one step of the basic 
circuit shown in Fig. Pfa) and S is the von Neumann entropy 
of the system. 
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FIG. 6: Plots of ST*, S, and S° ut -S for 4 qubits as functions 
of e. ST is the effective entropy of them after one step of 
the circuit shown in Fig. 

that of the 3-qubit basic circuits. It, however, induces 
a larger amount of correlation in comparison with the 
3-qubit basic circuits. The discrepancy between the ef- 
fective entropy of the output, S° ut , and the von Neumann 
entropy S is shown in Fig. 

When one compares the mean discrepancy (S° ut — 
S)/rib (rib — 3,4 is the number of wires) for the 3-qubit 
circuit with that for the 4-qubit circuit, as shown in Fig. 
it is clear that the increase in the amount of correla- 
tion in one time step is so large that the 4-qubit circuit 
should not be used as a basic boosting operation. This 
is the reason why we have chosen the 3-qubit circuits. 



V. AUTOMATIC CIRCUIT GENERATION 

A molecular simulation is applied for designing a quan- 
tum initialization circuit in the boosting scheme on con- 
ventional computers. It is implemented as a circuit com- 
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Original bias e 

FIG. 7: Plots of (S° ut " S)/n b for the 3-qubit circuit (n b = 3) 
and the 4-qubit circuit (nj, = 4) as functions of the original 
bias e. 



GEN_INIT_CIRCUIT 

CREATE_SYSTEM 
* 

Rearrange the pointers to the qubits. 

BOOST in every three pointers: 

Apply the basic boosting 
circuits to the three pointers. 



UPDATE_BIAS 

Output the circuits iff the 
boosting is successful; 
otherwise, undo the circuits 
and the biases. 



L 



FIG. 8: The brief flowing chart of the simulation algorithm 
of Program 1 to generate an initialization circuit on the basis 
of the boosting scheme. We show the program in Appendix 
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poser for the scheme, offering the forecasts of biases dur- 
ing the initialization. The details of the algorithm for 
the designing are described in Program 1 of Appendix 
IO In this section, instead of the long program, the brief 
flowing chart of it is shown in Fig. [3] It utilizes a virtual 
molecular system to mimic an n-qubit ensemble as illus- 
trated in Fig. ^ As we continue to use the same basis 
{|0 ■ • • 0), . . . , |1 • ■ • 1)}, superpositions of multiple compu- 
tational basis vectors are not produced by the boosting 
operations. Hence a molecule is represented as an ar- 
ray of n binary digits in the simulation. We consider N 
molecules that consist of n spins with bias £j in the ith 
spin (N > n). 

The outline of Program 1 is as follows: A table of point- 
ers to qubits is used instead of direct handling of them. 



We consider the subgroups which consist of every three 
pointers from the upper side to the lower side in the table, 
although those are not explicitly written in the program. 
Let frandQ be a function which returns a pseudo ran- 
dom real number between and 1. This is called n-N 
times in the procedure CREATEJ3YSTEM to set up 
the virtual molecular system. At each call of frand(), the 
bit representing the ith spin in the mth molecule is set 
to when frandQ < (1 +£j(0))/2 and set to 1 otherwise. 
The desired original biases £i(0), . . . , e„(0) are prepared 
in this way. Then the pointers are arranged in the or- 
der of the biases, from the larger to the smaller, in the 
table (from the upper side to the lower side). The main 
part of this simulation after the setup is an iteration in 
the procedure GEN JNIT.CIRCUIT. This iteration 
contains the following operations: Firstly, BOOST is 
applied to every subgroup mentioned above. BOOST is 
the boosting operation for three qubits a, b, and c. We 
apply the circuit of Fig. [2Ja) to them. When the bias 
of b becomes negative, the circuit of Fig. $[b) is also ap- 
plied. The boosting circuits applied to the three qubits 
are written out if the bias of a is increased by them; 
otherwise, we undo the boosting. The biases after the 
boosting are calculated by using UPDATE_BIAS in 
each call of BOOST. Secondly, we rearrange the table 
so that the pointers to qubits with larger biases appear in 
the upper side of it and those with smaller in the lower 
side. The iteration of above operations is terminated 
when cold qubits are not newly produced in the last s t 
steps or when the depth of generated circuits reaches a 
specified terminal depth of circuits. After this iteration, 
qubits with very large biases are pointed from pointers 
located near the top of the table; the component state 
|0)® with a very large population is achieved by collect- 
ing them. 

This algorithm assumes that a CNOT operation can 
be used with any two qubits in the system. When this 
is not guaranteed for dispersed qubits, SWAP operations 
are used to gather them into one location in the molecule. 
This change is easily realized by arranging the qubits 
with SWAP gates in GEN_INIT_CIRCUIT instead of 
arranging only the pointers. 

The complexity of the algorithm is 0(n • N) in space 
and 0(n 2 ■ N) in time in conventional computers. For 
a practical use of the simulation for n < 10 4 , it would 
be proper to set N ~ 10 4 n and St = 5 + [O.lnJ. Our 
simulation can be used for generating the circuits up to 
10 4 qubits for arbitrary biases within 1 terabit of memory 
space. 

In addition, the present algorithm is equivalent to the 
original one proposed by Schulman and Vazirani in Ref. 
0. In the iteration of the original algorithm, at each 
step, qubits are rearranged so that qubits with large bi- 
ases appear first followed by qubits with relatively small 
biases, discarding those with very low biases. Note that 
we do not need to discard qubits with very low biases 
actually. They are dispersed far away from those with 
large biases by the rearrangement automatically. In the 



present algorithm, in the iteration, the pointers to qubits 
are rearranged in the order of the biases. Therefore, it is 
equivalent to the original one. 



Example A 

An example ^(| of the output circuit from the program 
of Program 1 is shown in Table ITT1 It is generated for the 
7-pin system, 1 H- 1 H- 1 H- 1 H- 1 H- 1 H- 1 H with the uniform 
bias 0.6. Figure [U] shows that, according to the numeri- 
cal matrix calculation result, the bias of the first spin is 
enhanced as it is designed to. Although it appears to be 
possible to enhance the bias further at a glance, it is im- 
possible for the boosting scheme owing to strong classical 
correlations among the qubits. An NMR experiment of 
this bias boosting will be hopefully possible in future. 



TABLE II: The initialization circuit generated for the 
spin system, 1 H- 1 H- 1 H- 1 H- 1 H- 1 H- 1 H. The spins are labeled 
1, . . . , 7, from left to right, respectively. The transformations 
in this table are executed from left to right and upper to 
lower. X, CN, and Fr read NOT gate, controlled-NOT gate, 
and Fredkin gate (controlled-SWAP gate), respectively, fol- 
lowed by numbers in parentheses indicating which qubits are 
handled. In CN(a,fo), the control bit is a and the target bit is 
b; in Fr(a b, c), the control bit is c and the target bits are a 
and b. 



CN(2,3);X(3);Fr(l 2, 3) 

CN(5,6);X(6);Fr(4 5, 6) 

CN(2,6);X(6);Fr(7 2, 6) 

CN(2,3);X(3);Fr(6 2, 3) 

CN(6,5);X(5);Fr(3 6, 5) 

CN(5,3);X(3);Fr(6 5, 3) 

CN(4,7);X(7);Fr(l 4, 7) 



X(3) 
X(6) 
X(6) 
X(3) 
X(5) 
X(3) 
X(7) 



1 H 1 H 1 H 1 H 1 H 1 H 1 H 



1 H 1 H 1 H 1 H 1 H 1 H 1 H 




(b) 



FIG. 9: The result from a matrix calculation for the first 
example, (a) The original biases of the 7-spin system, (b) 
Redistributed biases after the operations of Table [H] One of 
the J H spin biases is boosted under the preset low tempera- 
ture such that the original bias of the hydrogen is 0.6, in this 
matrix calculation. 



Example B 

The second example is an output circuit for qubits 
with initially nonuniform biases. The circuit shown 
in Table IIIII is generated for the 9-spin system, 
iH-iH-iH-iH-^C-^C-^C- 31 ?- 31 ?. We used a modified 
program to boost mainly the bias of the 5th spin. The 
numerical matrix calculation of the boosting with the cir- 
cuit results in an increase in the bias of one of the 13 C 
spins as expected from the molecular simulation (Fig. 
Ill) [I . Of course, the effect of the boosting is too small to 
obtain an initialized qubit in such a small molecule. 



TABLE III: The initialization circuit generated for the spin 
system, 'H-'H^H-^-^C- 1 ^- 1 ^- 3 ^- 31 ?. The spins are la- 
beled 1, ... ,9, from left to right, respectively. The notation 
of the operations in this table is the same as that in Table ITT1 



CN(6,7) 
CN(9,2) 
CN(6,9) 
CN(1,4) 
CN(6,1) 
CN(4,3) 
CN(4,3) 
CN(7,8) 



X(7) 
X(2) 
X(9) 
X(4) 
X(l) 
X(3) 
X(3) 
X(8) 



Fr(5 6, 
Fr(8 9, 
Fr(7 6, 
Fr(2 1, 
Fr(9 6, 
Fr(7 4, 
Fr(9 4, 
Fr(5 7, 



7) ;X(7) 

2) ;X(2) 
9);X(9) 
4);X(4) 

i);X(i) 

3) ;X(3) 
3);X(3) 

8) ;X(8) 



X(4); 
















h 




T 


1 



(b) 



FIG. 10: The result from a matrix calculation for the second 
example, (a) The original biases of the 9-spin system, (b) 
Redistributed biases after the operations of Table ITTT1 One of 
the 13 C spin biases is boosted under the preset low tempera- 
ture such that the original bias of the hydrogen is 0.8, in this 
matrix calculation. 



VI. RESULTS FROM SIMULATION 

The simulation of Program 1 has been conducted under 
several different conditions of the original biases and the 
number of qubits, with the number of virtual molecules 
N = 5.0 x 10 6 . This value of N is large enough to ensure 
the reliability of the simulation as we show in Appendix 
lAl e co id was set to 2 x OJ 1 ^™ -5 ! — 1 (this is larger than 
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0.99 when n— S > 22) in Program 1. (S is the von Neu- 
mann entropy of a simulated system.) As described in 
Sec. Mil we propose to compute the sum of the binary 
entropies of individual qubits during the simulation be- 
cause it is a lower bound for the length of the compressed 
string at a particular time step t, at least approximately, 
in the boosting scheme. Recall that the effective entropy 
is thus the total binary entropy given by 



S e (t) 



tc- 



i=l 



log 2 



1 



log 2 



(e, = Si(t)) and that S e can be increased by the growth 
of classical correlations among qubits in the scheme. At 
the last step in the simulation, we ignored the preset 
value of £ co id and picked up / cold qubits so that the 
condition of ProbjOi • • • 0;} > 0.9 was satisfied; although 
this does not affect the initialization process itself, nor 
does it affect the time evolution of the effective entropy. 
Here we adopt the circuit of Fig.^a) as the unit quantum 
circuit for measuring the depth of circuits constructed by 
the scheme. In addition, let e denote uniform original 
bias for qubits. We assume that < e < 1 because a 
negative original bias can be inverted in advance. 

The logs recorded during each run of the simulation 
show that although the von Neumann entropy S — 
— Tr(/olog 2 p) is preserved, the number of qubits initial- 
ized after the boosting operations is much less than n — S 
because S e is increased in the early steps of the boost- 
ing procedure. This is well illustrated as the plots of 
S e against the depth of the circuits (denoted by d) for 
1000 qubits with originally uniform biases, as shown in 
Fig. ^2 F° r an example, let us take e — 0.7 in the fig- 
ure. Although the value of the von Neumann entropy is 
609.8, the value of S e at d = 40 is 806.8. The number of 
qubits that can be initialized is reduced from 390 to 193 
approximately. 

The relation of S* nd (the terminal value of S e ) to S 
was found to be of interest. We plot S^ nd /n against S/n 
for several different values of n in Fig. ^] The original 
biases were set to be uniform and the data were taken 
for < £ < 0.975. It suggests 



s: 



end . 



(31) 



for n > 70. In addition to this relation, it suggests that 
for n > 70 and < e < 0.65, 



Sl nd /n > 



(32) 



in consequence of the fact that a value of S^ nd /n was 
found to be insensitive to a value of n > 70 for this 
range of e. The figure also shows that the data points 
of Sl nd /n for n — 1000 is slightly closer to the curve of 
U S/n than those for smaller values of n. Furthermore, 
more than 60% of the increase in S e occurs during d < 5 
when n — 1000; the number of qubits to which a qubit is 
classically correlated is nearly 3 5 = 243 < 1000 for this 
depth. Hence the growth in S e /n during d < 5 is almost 
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FIG. 11: Plots of the effective entropy S e against the depth of 
circuits from the simulation logs for n = 1000. The simulation 
started with qubits with uniform biases. The data points were 
obtained for the different values of original bias: e = 0.3, 
£ = 0.5, e = 0.7, and e = 0.9. 



unchanged for larger values of n. Thus, as n becomes 
large, S e e nd /n varies only slightly, getting closer to \/Sjn. 
Therefore, Eq. (J3TJ for < e < 0.975 and inequality Q 
for < e < 0.65 are true for the circuits with larger 
widths, i.e., they are true for v n > 70. Beca use t he 
data points of Sf ld /n are in the interval of [yj S/n — 
0.044, y/Sjn~+ 0.032] when n = 1000, the relation of S c e nd 
to S for n > 1000 and < e < 0.975 is given by 



oend 



(33) 



where S is a real number (—0.05 < S < 0.04). This 
phenomenon is originated from the fact that the basic 
circuits increase the effective entropy of 3 qubits input 
with independent uniform biases unless their value is 
or ±1 as we have already seen in Fig. In addition, 
the value of Sf Ld /n is not sensitive to the preset value of 
e co id as long as e co u > 0.92 for e < 0.5 and e co u > 0.99 
for e ^ 0.5 as shown in Fig. El Therefore, the above 
relations, inequality Il-Sl'l) and Eq. (|33|l . are unchanged 
even if one presets e co id in another way at one's option 
as long as e co id is suitably large. 

Moreover, the initialization rate l/n is not rapidly im- 
proved by increasing the value of n as shown in Fig. 1141 
where / is the number of cold qubits picked up at the 
end of the program. Although / is expected to be close 
ton - Sl nd for large values of n, only I ~ Q.7(n — S^ nd ) 
is achieved even when n — 1000 and e = 0.9. This is 
due to the increase in the threshold of the bias for cold 
qubits picked up at the terminal step, since it becomes 
large for large n in order to avoid reduction in the pop- 
ulation of the component state |0i • • • 0;). There are still 
some other plots of data: The maximal efficiency of ini- 
tialization, r e = (n - Sl nd ) I '(n - S) , is plotted in Fig. H31 
It is not large enough unless the original bias is almost 1 
if we regard the scheme as a classical data compression 
method. Matrix calculations were conducted to verify 



10 




FIG. 12: Plots of St nd /n against S/n for the different values 
of n: n = 70, n = 200, and n = 1000. The values of original 
uniform bias e are also shown on the horizontal axis. Each 
run of the simulation was terminated at the depth of d — 100. 
The plots are data points from the simulati on, all of which 
approximately lie on the curve of Sl" d /n — \J S/n. The facing 
arrow pairs indicate the most distant data points from the 
curve in the vertical direction for n = 1000, accompanied by 
the distance values (both upper and lower side). 
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FIG. 14: The initialization rate l/n against n for the values of 
uniform original bias: e = 0.3, e = 0.5, e = 0.7, and e = 0.9. 
I is the number of cold qubits picked up at the last step in 
the simulation. For each data point, we chose the largest of 5 
samples of / that were calculated with different random seeds. 
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FIG. 13: Plots of S% /n against e co id for several different 
values of original uniform bias e. n was set to 100 for the 
data points in this figure. 
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FIG. 15: The maximal efficiency of initialization, r e = (n 
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curve is the spline fit 



en )/(n — S), for n = 1000 against the uniform original bias 
The crosses are data points from the simulation and the 



same amount of increase in Sl nd as that for originally 
uniform biases (see Appendix [BJ). 

One can conclude from the results of our molecular 
simulation, especially from Eq. (|33fl . that the number 
of initialized qubits after the application of the boosting 
scheme is 



(34) 



the behavior of the automatically generated circuits for 
n < 13, indicating the reliability of our molecular sim- 
ulation. Two of them have been already shown in the 
examples of Sec. [V] The reliability of the simulation is 
also statistically evaluated in Appendix El The relation 
of S e J ld to S was obtained also for the biases that are 
originally nonuniform when n — 70, showing almost the 



for n > 1000 and < e < 0.975, where S is the von 
Neumann entropy of the system, S e (—0.05,0.04), and 
(3 is a small number that has appeared in inequality i|23|) 
of Sec. 1 1 1 1 1 (3 is negligible for large I, in particular for 
n > 1000 and e ^ 0.3. When it is required that the bias 
of an initialized qubit picked up at the terminal step tends 
to 1 with increasing n, (3 = o(l). In addition, combining 
one of the results inequality l-'il'l) with inequality l|23|) , we 
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find that, for n > 70 and < e < 0.65, 

I < n (l - y/S/ri + 1 ). (35) 

VII. DISCUSSION 

Much work has been done on the theory of q uantum 
data compression of mixed states 0, 0, EE l26l I2H EsL 
Efll ]. Nevertheless, our numerical analysis is the first eval- 
uation of large quantum circuits actually generated by 
a compression scheme for ensemble quantum computing 
as far as the authors know. Calculating an actual com- 
pression rate for a specific case is quite another subject 
than a theoretical proof of a limit to a compression rate. 
For an initialization of ensemble computers, the string to 
be compressed is a macroscopic string which is a prob- 
abilistic mixture of microscopic strings. It is probably 
impossible to generate a quantum circuit providing the 
optimal compression rate for this case as long as only the 
biases of qubits are referred to at each step of the circuit 
designing. It would be possible to construct the optimal 
compression circuit when it is generated by referring to 
the populations of component states. So far, efficient and 
effective initialization schemes are those which target the 
typical sequences EE HE • 

The general difficulty in the evaluation of a quan- 
tum data compression scheme is caused by correlations 
among qubits after multiple-qubit operations. Our pro- 
gram has demonstrated that it does not cost so much 
resources to simulate a scheme for initial thermal states 
in which only permutations of component states are used 
since only classical correlations can be induced. The 
boosting scheme is one of this special type of compres- 
sion methods although its initialization rate is not large. 
There are many polynomial-time classical data compres- 
sion schemes offering the average code length approx- 
imately equal to the Shannon entropy of the original 
probability distribution for suitably large n, e.g., enu- 
merative coding |3l| . some of which were reconstructed 
as quantum algorithms [2E HE 13 HE • A common effi- 
cient classical data compression may be used for a more 
effective initialization as long as the number of required 
clean qubits of workspace is log 2 [0(poly(n))] [2E. With 
a proper molecular simulation, a quantum circuit con- 
structed by a classical data compression scheme would 
be easily evaluated. 

To realize a high initialization rate, the value of the 
total binary entropy S e at the terminal step in an initial- 
ization scheme must be nearly equal to that of the von 
Neumann entropy of the system. In the boosting scheme, 
as illustrated in Fig. ^2 although the growth of S e is a lit- 
tle suppressed, the suppression is not sufficiently strong. 
The amount of correlation S e — S grows step by step, 
mainly in the early steps, in the bias boosting process. In 
other words, a classical string whose macroscopic bits are 
strongly correlated to each other can only be compressed 
by the scheme without an increase in the amount of cor- 



relation, although such a string does not exist initially 
in a common NMR system. An effective initialization 
scheme must have an ability to suppress the growth in 
the total correlation. 

Finally, we have to discuss the curious relation of S^ nd 
to S that has been suggested by the simulation results. 
The relation given by Eq. I|31[) is simple despite the com- 
plicated process of the growth in S e . It is difficult to 
predict the growth without the molecular simulation. In- 
deed, there are nothing but a few trivial things that sup- 
port an intuitive interpretation of the relation: 

1. Sf ld = S when e = or e = 1. 

2. Sl nd increases monotonically as S increases, at least 
approximately. 

But there might be also a profound physics in the simple 
relation, although it is, of course, a result for a special 
scheme. We have calculated the growth of S e only in the 
circuits that comprise the basic circuits shown in Fig. |2 
It is of interest to examine a growth in the sum of binary 
entropies in a similar scheme with another basic circuits 
as well as that in other compression schemes. 

VIII. CONCLUSION 

A molecular simulation has been demonstrated to eval- 
uate the effectiveness of Schulman and Vazirani's boost- 
ing scheme for the initialization of NMR quantum com- 
puters. We have confirmed that a generated quantum 
circuit will enhance the biases of target qubits as it is de- 
signed to. However, our results have also indicated that 
the scheme is rather inefficient with respect to the ini- 
tialization rate even for a large number of qubits. When 
we use 3-qubit basic circuits, the rate is at most ap- 
proximately 1 — \/Sjn even for n > 10 3 as long as 
< £ < 0.975, where S is the von Neumann entropy 
of a whole n-qubit ensemble. This is owing to a large 
increase in the sum of the binary entropies of individual 
qubits. This increase means that a large amount of clas- 
sical correlation is induced in the macroscopic n-qubit 
string. It is to be hoped that an advanced algorithm 
that suppresses the growth in the amount of correlation 
will be constructed to improve the rate. 
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APPENDIX A: RELIABILITY OF THE 
SIMULATION 

Because our simulation uses a limited number of vir- 
tual molecules, there are, to some extent, errors in the 
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output data. The errors comprise both systematic and 
statistical errors, since the time evolution of the virtual 
molecular system is deterministic in Program 1 while the 
setup of the system is dependent on the random seed (see 
also Appendix [Ujl . Here, we present a statistical evalua- 
tion of the reliability of the simulation system although 
the errors are expected to be small in the data used in 
Sec. I VII as we have set the number of molecules suitably 
large. 

In order to evaluate the reliability, the average of sam- 
ple values of S^ nd was plotted against the number of 
molecules, N, as shown in Fig. Sf Ld is the value 
of the effective entropy at the depth of d = 100 as no 
further increase in S e was found at larger depth. Con- 
sidering the running time of the simulation, we chose the 
original uniform bias e — 0.5 and the number of qubits, 
n = 100. The data were obtained by executing the pro- 
gram 60 times for each value of N with different random 
seeds. In the figure, the error bars represent 99% con- 
fidence intervals associated with each mean. The half 
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FIG. 16: The average values of S^" d (solid squares) against 
the number of virtual molecules N when e = 0.5 and n = 100. 
The error bars represent 99% confidence intervals associated 
with each mean. 
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FIG. 17: Plots of V s against the number of virtual molecules 
N. V s is the variance of the sample values of St nd at each 
value of N. The sample values are the same as those used in 
Fig.EO 



APPENDIX B: INITIALIZATION FOR 
NONUNIFORM ORIGINAL BIASES 

We believe that the initialization for qubits with 
nonuniform biases is as effective as that for qubits with 
uniform biases, in the boosting scheme. Consider n 
qubits in the structure • • • -A-B-A-B-A-B- • • • where n 
is even; A and B represent the atom with original bias 
£a > and that with eb > 0, respectively. The initial- 
ization for this structure is equivalent to the one for the 
structure A- A- • • • -A-A-B-B- • • • -B-B since we arrange 
the pointers to the qubits in the order of the biases. 

When ea ^> Eb {i.e., < e_e <C 1), a pointer to A is 
rarely grouped with a pointer to B in early time steps in 
the scheme; the increase in the effective entropy S ej AB 
occurs separately in A- A- • • • - A- A and in B-B- • • • -B-B. 
Together with Eq. l|31fl . this leads to that the terminal 
value of S e .AB is given by 



oend 
°e,AB 



(Bl) 



width of the confidence interval is smaller than 0.4 bit 
for the values of N > 5.0 x 10 5 . Because the mean value 
of S^ nd was found to be insensitive to the value of N and 
we used the same number of samples for each data point, 
this indicates that the systematic errors are negligibly 
small for such large values of TV. Along with the growth 
of TV, the errors are reduced even further. Moreover, the 
variance of the samples, V s , is less than 0.1 bit 2 when 
N > 5.0 x 10 5 (Fig. U3). Each value of V s was calcu- 
lated with the 60 samples of at the corresponding 
value of N. This result shows that we do not need to 
employ a statistical average value of Sf nd , i.e., a single 
sample value of S^ nd is reliable enough, when N is set 
to such a large number. The non-averaged sample val- 
ues of Sl nd used in Sec. I VII are reliable since we have set 
N = 5.0 x 10 6 . 



where S A .. A < n/2 and S B . B < n/2 are the von Neu- 
mann entropies in the original distribution for A- ■ • • -A 
and B- • • -B, respectively. The increase in the effective 
entropy is smaller than or approximately equal to that 
for qubits with originally uniform biases, since 

oend / 



where S — S° A ri9 A 



norig 
B—B' 



On the other hand, when 
ea — £s, the increase in the effective entropy is almost 
the same as that calculated for the originally uniform 
biases. In this case, 



oend 
°e,AB 



(B3) 



Now we utilize the coefficient \ an d se t £ b = X s A 
to calculate the increase in the effective entropy in the 
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scheme when n — 70 for several different values of \ 
(0 < X < !)• The relation of S™jf B to 5 is well depicted 
in Fig. El The figure shows that, as long as 0.1 < x < 1, 
ST'ab i s n °t sensitive to the value of x an d all the points 
approximately lie on the curve of S^"f B = VnS. Hence, 
the use of Eq. I|B3|) is justified for n = 70. One can 
conclude that the number of initialized qubits is 



Iab 



< 



(l - \fS/n) 



(B4) 



for any realistic value of \ (0.1 < \ < 1). Although the 
relation has been verified in the structure that comprises 
two species of qubits, similar results are expected to be 
given for molecules that consist of atoms with many dif- 
ferent original biases. 
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FIG. 18: Plots of St n AB against S for X = 0.1, x = 0.4, 
X = 0.7, and x = 1 when n = 70. 



larger bias appears upper side and that with a 
smaller appears lower side. 

while new cold qubits are obtained in last St steps 
and depth of output circuit < preset max- 
imum depth; do 
j <— minimum j such that £tblh*l < Scold 
while j < n - 2 do 

call BOOST(tbl[j], tbl[j + 1], tbl[j + 2]) 
if £ tbi[j] is increased then 

Output the circuit applied in the last 
call of BOOST procedure with the 
indication of handled qubits, tbl[j], 
tbl[j + 1], and tbl[j + 2]. 

else 

Undo the last BOOST, 
endif 

3 «- 3 + 3 
endwhile 

Re-arrange tbl[ ]. 

Take a log of tbl[ ] and e tb i[ i for simulation use. 
endwhile 



end 



procedure CREATE.SYSTEM: 

Create N molecules. (A molecule consists of n binary 

digits.) 

begin 

for m <— 1 to N do 
for i <— 1 to n do 

if frand() < (1 + e 4 ) / 2 then 
ith bit of mth molecule <— 

else 

zth bit of mth molecule <— 1 
endif 
endfor 
endfor 

end 



APPENDIX C: DETAILS OF THE SIMULATION 
ALGORITHM 

Here, we show the details of the simulation algorithm 
for the bias boosting scheme in the following program of 
Program 1, which has been explained in Sec. IV! 

PROGRAM 1: The simulation algorithm for generating 
an initialization circuit based on the boosting scheme. 
Here, frandQ is a function that returns a pseudo random 
real number between and 1. A qubit with a bias 
greater than e co id is regarded as a cold qubit. 

procedure GEN_INIT_CIRCUIT: 
begin 

call CREATE_SYSTEM 

Create tbl[ ], a table of pointers to individual qubits 
(tbl[l] is the uppermost end). 

Arrange tbl[ ] so that the pointer to a qubit with a 



procedure BOOST(a, b, c): 
begin 

Consider the trio, ath bit, 6th bit, and cth bit. 
for m <— 1 to N do 

Apply Fig. [5{a) to the trio on mth molecule, 
endfor 

call UPDATE_BIAS(a) 
call UPDATE_BIAS(&) 
call UPDATE_BIAS(c) 
if £b < then 

for m <— 1 to N do 

Apply Fig. Hfb) to the trio on mth 
molecule, 
endfor 

call UPDATE_BIAS(&) 
endif 

end 

procedure UPDATE_BIAS(i): 
begin 
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s <— endfor 
for m <— 1 to N do p <— s j N 

if ith bit of mth molecule is then e\ *— 2 * p ■ 

s •*— s + 1 end 

endif 
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